Circulating Current States in Bilayer Fermionic and Bosonic Systems 
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It is shown that fermionic polar molecules or atoms in a bilayer optical lattice can undergo the transition to 
a state with circulating currents, which spontaneously breaks the time reversal symmetry. Estimates of relevant 
temperature scales are given and experimental signatures of the circulating current phase are identified. Related 
phenomena in bosonic and spin systems with ring exchange are discussed. 
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Introduction.- The technique of ultracold gases loaded 
into optical lattices fHO] allows a direct experimental study of 
paradigmatic models of strongly correlated systems. The pos- 
sibility of unprecedented control over the model parameters 
has opened wide perspectives for the study of quantum phase 
transitions. Detection of the Mott insulator to superfluid tran- 
sition in bosonic atomic gases , of superfluidity 
and Fermi liquid JUl in cold Fermigases, realization of Fermi 
systems with low dimensionality ]9l |l0[| mark some of the re- 
cent achievements in this rapidly developing field [11]. While 
the atomic interactions can be treated as contact ones for most 
purposes, polar molecules lfl2i [Til [l4Tl could provide further 
opportunities of controlling longer-range interactions. 

In this Letter, I propose several models on a bilayer op- 
tical lattice which exhibit a phase transition into an exotic 
circulating current state with spontaneously broken time re- 
versal symmetry. Those states are closely related to the "or- 
bital antiferromagnetic states" proposed first by Halperin and 
Rice nearly 40 years ago [15], rediscovered two decades later 
171 11811 and recently found in numerical studies in ex- 
tended t-J model on a ladder [19] and on a two-dimensional 
bilayer ll20ll . Our goal is to show how such states can be real- 
ized and detected in a relatively simple optical lattice setup. 

Model of fermions on a bilayer optical lattice.- Consider 
spin-polarized fermions in a bilayer optical lattice shown in 
Fig-Q] The system is described by the Hamiltonian 

H = V^2ni jr n2, r + ^2 ^2 V^ al n a ,rn a ',r' (1) 

r era' {rr') 

r a (rr') 

where r labels the vertical dimers arranged in a two- 
dimensional (2d) square lattice, er = 1, 2 labels two layers, 
and (rr') denotes a sum over nearest neighbors. Amplitudes t 
and t' describe hopping between the layers and within a layer, 
respectively. A strong "on-dimer" nearest-neighbor repulsion 
V ^> t, t' > is assumed, and there is an interaction between 
the nearest-neighbor dimers V aa , which can be of either sign. 

This seemingly exotic setup can be realized by using po- 
lar molecules IU3[[l4Tl . or atoms with a large dipolar magnetic 
moment such as 53 Cr lfl2ll . and adjusting the direction of the 
dipoles with respect to the bilayer plane. Let 9, ip be the polar 
and azimuthal angles of the dipolar moment (the coordinate 



axes are along the basis vectors of the lattice, z axis is perpen- 
dicular to the bilayer plane). Setting cp 



= ±- 



± 



:',jv 



ensures 



the dipole-dipole interaction is the same along the x and y di- 
rections. The nearest neighbor interaction parameters in (Q]l 
take the following values: V — (d^/l^)(l — 3 cos 2 9), and 
V{ 2 = V21 = (do/fl 3 ){l-3i2 -2 (4| cos (9+^ sin (9 cos v?) 2 }, 
Vd = V22 — V{ 2 = 0), where do is the dipole moment 
of the particle, £± and £\\ are the lattice spacings in the direc- 
tions perpendicular and parallel to the layers, respectively, and 
R 2 = £ 2 + l\ . The strength and the sign of interactions V, 

V' can be controlled by tuning the angles 9, ip and the lattice 
constants £j_, in. Below we will see that the physics of the 



problem depends on the difference 

v' = v[ x - v{ 2 , 



(2) 



with the most interesting regime corresponding to V < 0. 

Consider the model at half-filling. Since V ^> t, t', we may 
restrict ourselves to the reduced Hilbert space containing only 
states with one fermion per dimer. Two states of each dimer 
can be identified with pseudospin--i states |f) and ||). Second- 
order perturbation theory in t' yields the effective Hamiltonian 



Us = Y,{ J ( S r S r'+S V r S V r ,) + J z S^,}-HY,Si. 
(rr') r 



J = 4(t') /V, J Z = JA = J + V', H = 2t, 



(3) 



describing a 2d anisotropic Heisenberg antiferromagnet in 
a magnetic field perpendicular to the anisotropy axis. The 
twofold degenerate ground state has the Neel antiferromag- 
netic (AF) order transverse to the field, with spins canted to- 
wards the field direction. The AF order is along the y axis for 
A < 1 (i.e., V < 0), and along the z axis for A > 1 (V > 0). 
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FIG. 1: Bilayer lattice model described by the Hamiltonian QJ. The 
arrows denote particle flow in the circulating current phase. 



2 



The angle a between the spins and the field is classically given 
by cos a = H/ (2ZJS), where S is the spin value and Z = 4 
is the lattice coordination number. This classical ground state 
is exact at the special point H = 2SJy / 2(I + A) $2^. The 
transversal AF order vanishes above a certain critical field H c ; 
classically H c = 2ZJS, and the same result follows from the 
spin-wave analysis of (01 (one starts with the fully polarized 
spin state at large H and looks when the magnon gap van- 
ishes). This expression becomes exact at the isotropic point 
A = 1 and is a good approximation for A close to 1. 

The long-range AF order along the y direction translates 
in the original fermionic language into the staggered arrange- 
ment of currents flowing from one layer to the other: 



N y = (-Y(sy)^(-r(--(al 



-,f 



*i,r)>- ( 4 ) 



In terms of the original model ([T), the condition H < H c for 
the existence of such a staggered current order becomes 



t < 8(t r ) 2 /V. 



(5) 



The continuity equation for the current and the lattice symme- 
try dictate the current pattern shown in Fig.Q] This circulating 
current (CC) state has a spontaneously broken time reversal 
symmetry, and is realized only for attractive inter-dimer inter- 
action V' < (i.e., the easy-plane anisotropy A < 1) 12211 . 
If A = 1, the direction of the AF order in the xy plane is 
arbitrary, so there is no long-range order at any finite temper- 
ature. For A > 1 (i.e., V' > 0) the AF order along the z axis 
corresponds to the density wave (DW) phase with in-layer oc- 
cupation numbers having a finite staggered component. 

The phase diagram in the temperature-anisotropy plane is 
sketched in Fig. [2] At the critical temperature T = T c the dis- 
crete Z2 symmetry gets spontaneously broken, so the corre- 
sponding thermal phase transition belongs to the 2d Ising uni- 
versality class (except the two lines A = 1 and H = where 
the symmetry is enlarged to U(l) and the transition becomes 
the Kosterlitz-Thouless one). Away from the phase bound- 
aries the critical temperature T c ~ J, but at the isotropic point 
A = 0, H = it vanishes due to divergent thermal fluctua- 
tions: for 1 — A <C 1 and H <C J, it can be estimated as 



T c ~ J/m[min(|l- Ap\ J 2 /ff 2 )]. 



(6) 




density wave 



FIG. 2: Schematic phase diagram of the model Q}, ([3} at fixed H = 
It. The line A = 1 corresponds to the Kosterlitz-Thouless phase, 
with the transition temperature Tkt oc J / \n(J/H) at small H. 
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FIG. 3: Noise correlation function Q(r, r') from time-of-flight im- 
ages in the circulating current (CC) phase, shown as the function of 
the relative distance Q(r) — Q(r'), with Q(r) — mr/(ht) ex- 
pressed in units: (a) Q(r) = (0,0); (b) Q(r) = f (1,1). 
Changing the initial point Q(r) leads to the change of relative weight 
of the two systems of dips, which is the fingerprint of the CC phase. 



The quantum phase transition at T = 0, H = H c is of the 3d 
Ising type (except at the U(l)-symmetric point A = 1 where 
the universality class is that of the 2d dilute Bose gas rf23ll ). so 
in its vicinity the CC order parameter N y oc (H c — H)" with 
(3 ~ 0.313 [24], and T c oc JN? oc J(H C - H) 2f} . At T > T c 
or H > H c the only order parameter is (S x ), corresponding 
to the Mott phase with one particle per dimer. 

Bilayer lattice design and hierarchy of scales.- The bi- 
layer can be realized, e.g., by employing three pairs of mutu- 
ally perpendicular counter-propagating laser beams with the 
same polarization and adding another pair of beams with 
an orthogonal polarization and additional phase shift S, so 
that the resulting field intensity has the form [E± (cos kx + 

cosky)+E z coskz] 2 +E 2 cos 2 (kz+S). Taking^ = f (1+0 
and E% > E z (2E± +(E Z ), with £ = ±1 for blue and red de- 
tuning, respectively, one obtains a three-dimensional stack of 
bilayers, separated by large potential barriers U^d- Eq. <(3j im- 
plies V 3> t' ^> t, \ V'\, which can be achieved by making the 
z-direction potential barrier U± inside the bilayer sufficiently 
larger than the in-plane barrier Un , so that the condition t <C t' 
will be met; e.g., E z /Ej_ w 20, E z /E^ w 15 yields the bar- 
rier ratio Usd '■ U± : Uh of approximately 16 : 8 : 1, and the 
lattice constants £± w 0.45 A, t\\ = A, where X = 2ir/k is the 

laser wave length. The parameter V' has a zero as a function 
of the angle 8, so it can be made as small as needed. Taking 
A = 400 nm, one obtains an estimate of T c = (0.1 -j- 0.3) /j,K 
for cyanide molecules C1CN and HCN with the dipolar mo- 
ment do « 3 Debye, while the Fermi temperature for the same 
parameters is Tf rj (0.6-fl.3) /iK. This estimate corresponds 
to the maximum value of T c ~ J reached when V 1 ~ — J 
and t < J. The hopping t' was estimated assuming the in- 
plane potential barrier Un is roughly equal to the recoil energy 
E r = (hk) 2 /2m, where m is the particle mass. 

Experimental signatures. - Signatures of the ordered 

l2r3l in time-of-flight experi- 



phases can be observed ||25 



ments by measuring the density noise correlator Q(r, r') 
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(n(r)n(r')) — (n(r))(n(r')). If the imaging axis is perpen- 
dicular to the bilayer, n(r) = J2<?( a t r a a r) i s me l° ca l net 
density of two layers. For large flight times t it is proportional 
to the momentum distribution riQ( r y, where Q(r) = mr/ht. 
In the Mott phase the response shows fermionic "Bragg dips" 
at reciprocal lattice vectors g = (2Trh/£n ) 2'Kk/£u), 



G M {r,r') oc /o(r,r') 



-2{S*) 2 yj{Q(r) 



g 



■Q(r') 



In the CC and DW phases the noise correlator contains an 
additional system of dips shifted by Qb = (n / t\\,~K / 

g C c,uw(r,r') oc / (r, r') - 2{(S Z ) 2 + (S y ) 2 
x [1 + ( cos (Q x (r)£ {[ ) + cos (Q y {r)l^) 2 } } 
x J2s(Q(r)-Q(r')-Q B -g) (7) 

9 

In the DW phase (S z ) ^ 0, (S y ) = 0, and so the density cor- 
relator depends only on r — r'. In the CC phase (S z ) = 0, 
(S v ) ^ 0, and the relative strength of the two systems of dips 
varies periodically when one changes the initial point r, see 
Fig. |3] This Q-dependent contribution stems from the intra- 
layer currents (a^ r a a r ,) = (-) <T (-) r 5( rr /)i(5 y }/4, where 
| comes from the fact that the inter-layer current splits into 
four equivalent intra-layer ones (see Fig. [TJ, S/ rr >\ means r 



(rr') 
= piQl 



de- 



and r' must be nearest neighbors, and (— ) r 
notes an oscillating factor. If the correlator is averaged over 
the particle positions, the CC and DW phases become indis- 
tinguishable. A direct way to observe the CC phase could be 
to use the laser-induced fluorescence spectroscopy to detect 
the Doppler line splitting proportional to the current. 

Bosonic models.- Consider the bosonic version of the 
model (Q]l, with the additional on-site repulsion U. The ef- 
fective Hamiltonian has the form (O with J = — A(t') 2 /V 
and J z = V' + A(t') 2 {l/V - l/U). Due to ferromag- 
netic (FM) transverse exchange, instead of spontaneous cur- 
rent one obtains the usual Mott phase. CC states can be in- 
duced by artificial gauge fields 02711 : The vector potential 
A(x) = j-(x + 1/2) makes hopping along the x axis imagi- 

nary, t' i— > it'. The unitary transformation Sp y \— * ( — ) r S*' y 
maps the system onto a set of FM chains along the x axis, 
AF-coupled in the y direction and subject to a staggered field 
H = 2t along the x axis in the easy (xy) plane. In the ground 
state net chain moments are arranged in a staggered way along 
the y axis, so a current pattern similar to that of Fig. [Demerges, 
now staggered along only one of the two in-plane directions. 

A different type of CC states, with orbital currents localized 
at lattice sites, can be achieved with p-band bosons Q28ll . 

Yet another way to create a CC state in a bosonic bilayer is 
to introduce the ring exchange on vertical plaquettes: 



In pseudospin language, the ring interaction modifies the 
transverse exchange, J i— > J + K , so for K > one can 
achieve the conditions J > 0, J > \J Z \ necessary for the CC 
phase to exist. However, engineering a sizeable ring exchange 
in bosonic systems is difficult (see [29] for recent proposals). 

Spin-^ bilayer with four-spin ring exchange.— Consider 
the Hubbard model for spinful fermions on a bilayer shown 
in Fig. [U with the on-site repulsion U and inter- and intra- 
layer hoppings t and t', respectively. At half filling (i.e., two 
fermions per dimer), one can effectively describe the system 
in terms of spin degrees of freedom represented by the opera- 
tors S = jo^a^apCLp- The leading term in t/U yields the AF 
Heisenberg model with the nearest-neighbor exchange con- 
stants Jj_ = At 2 /U (inter-layer) and J\\ = A(t') 2 /U (intra- 
layer), while the next term, with the interaction strength J4 ~ 
10t 2 (t') 2 /U 3 , corresponds to the ring exchange i30ll3lll : 

H 4 = 2J 4 ^{(Si-S 2 )(Si<-S 2 

□ 

+ (S x • Sv)(S 2 ■ Sv) - (Si • S2OOS2 • Sv)}, (9) 

where the sum is over vertical plaquettes only (the interaction 
for intra-layer plaquettes is of the order of (t') A /U 3 and is ne- 
glected), and the sites (1, 2, 2', 1') form a plaquette (traversed 
counter-clockwise). In the same order of the perturbation the- 
ory, the nearest-neighbor exchange constants get corrections, 

J± >-> Jr = J± + J4, J\\ Jl — J\\ + J4/2, 



Wring = \K (b\A,r'l> 
{rr') 



2,A, 



h.c). 



and the interaction Jd = \J± along the diagonals of verti- 
cal plaquettes is generated. Generalization for any 2d bipar- 
tite lattice built of vertically arranged dimers is trivial. Since 
J± ^> J|| , J4, we can treat the system as a set of weakly cou- 
pled spin dimers. The dynamics can be described with the 
help of the effective field theory B32 [] which is a continuum 
version of the bond boson approach 1 33] and is based on dimer 
coherent states \u,v) = (1 — u 2 — v 2 )\s) + J2j( u j + ^ v j)\tj)- 
Here |s) and \tj), j — (x,y,z) are the singlet and triplet 
states, and u, v are real vectors related to the staggered mag- 
netization {Si — S2) = 2u(l — u 2 — v 2 ) 1 / 2 and vector chi- 
rality (Si x S2) = v(l — u 2 — v 2 ) 1 / 2 of the dimer. Using the 
ansatz u(r) = (— ) r cp(r), v(r) = (— ) r x(r), passing to the 
continuum in the coherent states path integral, and retaining 
up to quartic terms in it, v, one obtains the Euclidean action 

A = J dTd 2 r{h( V -drX-X- d T <p) (10) 
+ (<^ 2 + X 2 )(Jr. - ZZJi/2) - Z[J^ 2 + J 4X 2 } 
+ {Z/2)[J\\{d k if) 2 + MdkX) 2 }+ZU 4 (cp, X )}, 
where the quartic potential U4 has the form 
Ua = (f 2 + X 2 )[J\\<P 2 + JiX 2 ]/2 

+ Mv 2 + x 2 ) 2 + (J\\ + J4)(vxx) 2 - (ID 

(8) Interdimer interactions Jm and J4 favor two competing types 
of order: while Jm tends to establish the AF order (tp ^ 0, 
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X = 0), strong ring exchange J 4 favors another solution with 
V = 0, x 7^ 0, describing the state with a staggered vector 
chirality. It wins over the AF one for J 4 > Jy, J 4 > ^ Jr, 
which for the square lattice (Z — 4) translates into 

J 4 > max(J||, Jj_/9). (12) 

On the line J4 = J|| the symmetry is enhanced from SU(2) 
to SU(2) x U(l), and the AF and chiral orders can coexist: a 
rotation {ip + ix) >— ► (f + ix) ela leaves the action invariant. 

The chiral state may be viewed as an analog of the cir- 
culating current state considered above: in terms of the 
original fermions of the Hubbard model, the z-component 
of the chirality (Si x S 2 )z = §{( a ix a 2i)( a 2T a iT) — 
(a\^a i:i )(a\^a 2 ^)} corresponds to the spin current (particles 
with up and down spins moving in opposite directions). 

Summary.- I have considered fermionic and bosonic 
models on a bilayer optical lattice which exhibit a phase tran- 
sition into a circulating current state with spontaneously bro- 
ken time reversal symmetry. The simplest of those models 
includes just nearest-neighbor interactions and hoppings, and 
can possibly be realized with the help of polar molecules. 
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